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We study cosmological constraints on the various accelerating models of the universe using the time evolution of the cosmological 
Tedshift of distant sources. The important characteristic of this test is that it directly probes the expansion history of the universe. 
In this work we analyze the various models of the universe which can explain the late time acceleration, within the framework of 
General Theory of Relativity (GR) (XCDM, scalar field potentials) and beyond GR {f{R) gravity model). 



1. Introduction 



The recent accelerated expansion of the universe is one of 
the most important discovery in the cosmology. Whether this 
^ observed acceleration is due to some new hypothetical energy 
O component with large negative pressure (dark energy) within 
O^the framework of General Theory of Relativity, or due to modi- 
;-H fication in the GR at the cosmological distances (modified grav- 
, ^^ ity), is not known. Therefore many models have been proposed 
|in literature to understand the origin and nature of this present 
CO acceleration |T]. 

Basically the cosmic acceleration aff'ects the expansion his- 
tory of the universe. Therefore to understand the true nature 
of this driving force, mapping of the cosmic expansion of the 
universe is very crucial [2]. Hence, we require various obser- 
vational probes in diff'erent redshift ranges to understand the 
expansion history of the universe. The observational tools for 
probing the cosmic acceleration broadly fall into two categories: 
Geometrical and Dynamical probes. A Geometrical probe deals 
with large scale distances and volume which include luminos- 
ity distance measurements of SNe la, angular diameter distance 
from first CMB acoustic peak, Baryon Acoustic Oscillations 
(BAO) etc.. A Dynamical probe investigates the growth of mat- 
ter density perturbations that give rise to the large scale struc- 
ture such as galaxies, clusters of galaxies etc. in the universe. 

Using supernovae as standard candles is a popular method 
of constraining the properties of dark energy. Though this method 
is very simple and useful in constraining the various dark en- 
ergy models, at present the luminosity distance measurements 
suff'ers from many systematical uncertainties like extinction by 
dust, gravitational tensing etc. On the other hand, measur- 
ing the expansion history from growth of matter perturbations 
also has its limitations. It requires prior information of exact 
value of matter density, initial conditions, cosmological model 
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etc. 10,0]. So the question arises, "Is there any probe which is 
simple, depends on fewer priors and assumptions?" The pos- 
sible probe is "Cosmological Redshift Drift" (CRD) test which 
maps the expansion history of the universe directly. 

The CRD test is based on very simple and straightforward 
physics. However, observationally it is a very challenging task 
and requires technological breakthroughs |5]. The most re- 
markable feature of this probe is that it measures the dynam- 
ics of the universe directly - the Hubble expansion factor. This 
property makes it very special and unique. Further, it assumes 
that the universe is homogeneous and isotropic at the cosmo- 
logical scales [for more details see ref. |6]]. The time drift 
of the cosmological redshift probes the universe in the redshift 
2 < z < 5, whereas the other cosmological tests based on SNe 
la, BAO, weak tensing, number counts of clusters etc. have not 
penetrated beyond z = 2. The other advantage of this tool is 
that it has controlled systematical uncertainties and evolution- 
ary eff'ects of the sources. 

The aim of this letter is to employ CRD test to constrain var- 
ious accelerating models both within the framework of GR and 
beyond GR. We have used simulated data points for redshift 
drift experiment generated by Monte Carlo simulations with the 
assumption of standard cosmological model (A CDM) as refer- 
ence iallJ]. We put constraints on the Quintessence models 
(based on scalar field potentials) like PNGB, inverse power law 
and exponential potentials. The dark energy parametrization 
which has a variable equation of state is also investigated. We 
also put a bound on the f{R) gravity models i.e., Starobinsky 
model using the Cosmological Redshift Drift test. 

The letter is organized as follows. In Section 2 we review 
the theoretical basis of Cosmological Redshift Drift test. We 
also present here methodology and data used for this work. Var- 
ious models which can explain late time accelerated expansion 
of the universe are described in Section 3. Last section contains 
a summary and discussion of results. 
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2. Cosmological Redshift Drift Test 

2.1. Theory 

A test which can trace the dynamical expansion history of 
the universe was proposed by Sandage |7|]. The expansion of 
the universe is expressed in terms of a scale factor, a{t). There- 
fore, the time evolution of the scale factor, or change in red- 
shift, z, directly measures the expansion rate of the universe. 
The redshift, z, of an object as determined today will be dif- 
ferent from its measured value after a time interval of several 
years. Sandage also stressed on the fact that the redshift drift 
signal, z, is very small. The significance of this tool has been 
discussed by several authors | 8]. Loeb was the first to sug- 
gest the possibility of measuring the redshift drift by observing 
Lyof absorption lines in the spectra of quasars (QSOs) |9]. This 
reinforced the importance of this probe. 

The observed redshift of a distant source is given by 



z{to) = — — - 1 
a(ts) 



(1) 



where ts is the time at which the source emitted the radiation 
and ^0 is the time at which the observation is made. In writ- 
ing the above expression, we ignore any peculiar motion of the 
object. The redshift of the source after the time interval of A^o 
becomes 

a(to -h A^o) 



z(to + A^o) 



1 



(2) 



a(ts + Ats) 

where Ats is the emission time interval for the source. In the 
first order approximation we can write 

Az_ _ (a(to) - a(ts)) 
Ato 

or 



a(ts) 



(3) 



H(z) 



Hn 



(4) 



The above equation is also known as McVittie equation 
This clearly shows that z traces H(z), which is the Hubble factor 
at redshift z. As stated earlier z measures the rate of expansion 
of the universe: z> and < indicates the accelerated and de- 
celerated expansion of the universe, respectively. For a coasting 
universe z = 0. The redshift variation is related to the apparent 
velocity shift of the source: 

Az 



Ay = c- 
Thus we can write 



cHq 



H(z) 



Hn 



(5) 



(6) 



where v = Av/A^o and Hq = 100 /z km/s/Mpc. In a standard 
cosmological model (ACDM), with a time interval of A^o = 10 
yr, the change in redshift is Az ~ 10"^, for a source at redshift 
z = 4. The corresponding shift in the velocity is of the order of 
Av 6 cm/s. To measure this weak signal, Loeb pointed out 
that observation of the Lya forest in the QSO spectrum for a 
decade might allow the detection of signal of such a tiny mag- 
nitude. 



2.2. Data 

In the near future, a new generation Extremely Large Tele- 
scope (ELT, 25 - 42 m diameter) equipped with a high reso- 
lution, extremely stable and ultra high precision spectrograph 
(CODEX) should be able to measure such a small cosmic sig- 
nal. The CODEX (COsmic Dynamics Experiment) operates 
in the spectral range of 400-680 nm with resolving power R 
= 150000. Several groups have performed Monte Carlo (MC) 
simulations of quasars absorption spectra and obtained 

the z measurements. In this work we have used three sets of 
data (8 points) for redshift drift experiments generated by MC 
simulations lldllSll. 

These three datasets are generated by three diff'erent ap- 
proaches 1 6]. The data set with error bars are generated by 
assuming the total duration of 20 years for observations and 
standard input cosmological model with Hq = 70 km/s/Mpc, 
flm = 0.3 and Da = 0.7. In every approach it is assumed that 
normalized observational set-up parameter, O, equal to 2. This 
parameter controls the telescope size, efficiency and integration 
time. In the first approach, the data points are selected by small- 
est value of CTy. In this approach set of 20 QSOs are distributed 
in the four equally sized redshift bin. The second set of data 
points are generated by selecting the targets by larger value of 
|v|/(Tv. In this case, Nqsos = 10 are distributed in two redshift 
bins. Finally the last two data points are generated by increas- 
ing the sensitivity of data towards Qa. 

In principle, redshift experiment involve the measurement 
of velocity shift between pair of spectra of same object sepa- 
rated by large time interval of many years. As mentioned ear- 
lier. There are many candidates available for this measurement 
but the most promising one is the absorption features of the Lya 
forest in the high redshift QSOs. The main reason is that the 
peculiar motion associated with intervening gas is negligible as 
compare to cosmic signal and numerous lines are present in the 
single spectrum. Using MC simulations, the desired statistical 
error on velocity shift can be written as 



cTy = 1.35 



S/N 
2370 



QSOs 

30 



-1/2 



1 + ZqsOs 



cm/s 



with a = 1.7 for z < 4 and a = 0.9 for z > 4, where S/N is 
the signal-to-noise ratio is of the order of 13500 per 0.0125 A° 
pixel, Nqsos is the total number of quasar spectra observed and 
Zqsos is their redshift. QSOs are the brightest sources which 
exist even at the redshift close to 6 and since we are observing 
it from ground, the existence of atmosphere put the lower limit 
of redshift z^2. Using full MC simulation to achieve a radial 
velocity accuracy of 3 cm/s, can be obtained with 3200 hours 
of observation with 42 m ELT. For more details see Ref. 

2.3. Method 

We perform the analysis to find the best fit values of the 
cosmological parameters and to find the bounds on them. 
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iyth{zk,Kp)-\Mc{zk)f 
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where p are the model parameters. Here i/th and vmc are the 
expected and the simulated values of the velocity drift respec- 
tively. The error bars on the velocity drift are denoted by 0-^. 
To draw the likelihood contours at 1, 2 and 3cr, Ax^ - ~ 
xl^i^ = 2.30, 6.17, 11.8 respectively in two-dimensional para- 
metric space. 

3. Models 

In this work, we investigate various models of the universe 
which can explain late time acceleration both within the frame- 
work of GR (Einstein gravity) and beyond GR (modified grav- 
ity). 

3.1. Models Based on Einstein Gravity 
3.1.1. XCDM Model 

One of the widely studied model of dark energy is XCDM 
parametrization. In this model the dark energy is characterized 
by an equation of state Wx = pip, where p and p are pressure 
and energy density respectively. Further, the equation of state 
does not evolve with time. For acceleration in this model Wx < 
-1/3. The Friedmann equation in this model is: 



H{z) 



(8) 



Here we assume the universe is flat i.e., Q.m + = 1- and 
Qx are the fractional matter and dark energy densities at the 
present epoch respectively. 

3.1.2. Scalar Field Cosmological Models 

In order to get Hubble parameter, //(z), in the scalar field 
cosmology, we have to solve the following equations of motion. 
The Einstein field equations can be expressed as: 



H 



(9) 



Here dots are derivative w.r.t. time. The scalar field equation of 
motion is 



= -3H4> - V\(p) 



(10) 



where prime stands for a derivative w.r.t. the scalar field 0, and 
the Hubble parameter is 



SttG 



Pm + Y + V((p) 



(11) 



Here is energy density of matter and nipi = G is the 
Planck mass. The equation of state oj is defined as 



OJ = 



02 - 2V((P) 



02 + 2V((f)) 

We analyse three scalar field potential models. 



(12) 



• Inverse power law potential: In this model the scalar field 



potential is of the form: 
k 



V(0) 



(13) 



where both k and y are positive and dimensionless con- 
stants. G is the gravitational constant 1IT4I] . This potential 
shows the tracking behaviour in which scalar field start 
from wide range of initial set of conditions and in the 
late time it approaches the cosmological constant. The 
Friedmann equation in terms of dimensionless parame- 
ters becomes [for details see ref. [15]] 



H 

Ho 



(14) 



where km^^JH^ - ?)6Klh?', with ^ > 0. The dimension- 
less variables are defined as: 



X : 



4^/7^ 



0, Y 



H 



Jnpi 



With these definitions we can write the equations of mo- 
tion in the following first order form 



dX 

dz 

dY 

dZ 
~di 



1 



(1+z) 
1 

1 



3X 



ykm^pi ^-(r+i) 



2Hl 



3 1 X^ kmli z-r 
2 ^^Y~MI^~ 



^0 



(15) 



(1+z) 



We solve the above system for parameter values k = Kx 
10-120^ ^ ^ 57 X 10-^^/z^y, mpi = 1.2 X lO^^eV, 
and h = 0.7. For these values for given set of parame- 
ters we choose initial values (redshift , z ~ 3000) Z/ = 
2.1, Xi = 10"^. Now for diff'erent values of the parameter 
K, we generate initial data sets by choosing ¥( in such a 
way that 7 = 1 at z = 0. 

The variation of velocity drift w.r.t. source redshift is 
shown in Fig. [T] 

Exponential potential: This potential has the following 
form 



y(0) = Voexp(-cp/f) 



(16) 



where Vq and / are positive constants The ex- 

ponential potential has the capability of producing scal- 
ing solutions which further scales the background energy 
density. But it require fine tuning of parameters to pro- 
duce the late time acceleration. The corresponding Hub- 
ble parameter in this model is 



(17) 
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Figure 1: Velocity drift (cm/sec/yr) versus redshift for inverse power law po- 
tential with different values of K (y is fixed at 0.001). The bold black curve 
correspond to the AC DM model. 



where j3 = m\l{^7if'^). The dimensionless variables are 



defined as: 



X = 0, Y : 



H 



t 
f 



In 



Vo 



Field equations can now be written as 



dX 
dz 
dY 
~d^ 
dZ 
dz 



1 



1 

1 



3X - STTyS— 



3^ 1 ^-z 

-7+ 47T — 

2 4/3 Y Y 
X 
Y 



(18) 



We solve the above system for parameter values Hq = 
1.67 X 10-^^/z eV, jfipi = 1.2 X 10^^ eV and h = 0.7. For 
these values for given set of parameters we choose initial 
values (redshift , z ~ 3000) Z, = 1.5, X, = 10"^ Now 
for diff'erent values of the parameter /3, we generate initial 
data sets by choosing 7/ in such a way that Y = 1 at z = 0. 

The variation of the velocity drift for exponential poten- 
tial model are shown in Fig. O It is clear from the figure 
that after certain value of j3, the universe remains in the 
decelerated phase. 

• PNGB Model : In this (Pseudo-Nambu-Goldstone Bosons) 
model the scalar field potential has the following func- 
tional form 



y(0) = (1 + cos((P/f)) 



(19) 



where M and / are positive constants lfl5[ [l7. [13] . The 
main motivation to study this potential is because of its 
special properties which not only explain late time ac- 
celeration but its ability to off'er solution to the cosmic 



Figure 2: Variation of velocity drift (cm/sec/yr) with redshift z for y = 
Voexpi-(f)/ f) . The bold black curve correspond to the AC DM model. 

coincidence problem. The Hubble parameter for PNGB 
model in the term of dimensionless parameters is 



A-TT 

Y^ = —a^ X^ + 1+ cos(Z) + a„ 



where 



X : 



H 



f 



(20) 



(21) 



and a = f/rupi. The full dynamical system in terms of 
variables X, Y and Z can be written analogues to earlier 
two cases as 



dX 
dz 
dY 
~di 
dZ 
~di 



1 



(1+z) 
1 

(1+^) 
1 



(1+z) 



mJ/sinZ 

3 oX^ 

X 
Y 



■ 4nm 



' 1 -h COSZ'^ 



(2: 



where = M^/(mpiHo)^. We choose initial values (red- 
shift , z ~ 3000) Zi = 1.5, Xi = 10-\ Now for diff'erent 
values of the parameter a, we generate initial data sets by 
choosing 7/ in such a way that Y = 1 at z = 0. 

The variation of the velocity drift for PNGB model are 
shown in Fig. [3l It is clear from the figure that after 
a certain critical value of a, universe always remain in 
the deceleration mode when the light is traveling toward 
us. For smaller values of a, it shows the universe has 
undergone the decelerating phase twice. 

3.2. Modified Gravity Model 

An important class of models which has attracted consider- 
able attention in past few years is the one that modifies Einstein- 
Hilbert action by replacing Ricci tensor by an arbitrary function 
of curvature 



m_ 

l6nG 



(23) 




Figure 3: Variation of velocity drift (cm/sec/yr) with redshift zforV = M^(l + 
cos(0//)) for M = 0.004 ^^heV. The black curve correspond to the AC DM 
model. 



In this letter we study the f(R) theory model given recently by 
Starobinsky, Hu and Sawicki ||20i|2l|] (see also, ||23] ) 



f(R) = ARo 



1 + 



1 



(24) 



where A, n and Rq are positive parameters. For necessary steps 
leading to calculation of Hubble parameter see Dev et al. [23]. 
These models can evade local gravity constraints and have the 
capability being distinguished from the cosmological constant. 

The Fig. [4] displays the variation of velocity drift for f{R) 
model. There is very small difference between the variation of 
H{z) with z for n =1 and n = 2 at the redshift z < 1.8 B . 
Hence for redshift z > 1.8, the variation of velocity drift be- 
comes independent of the value of n. 




observational tests either belong to distance based mtWxod^ such 
as SNe la luminosity distances, angular size of compact radio 
sources, BAO, CMBR, gravitational lensing etc. or time based 
methods like Absolute age method. Lookback time method and 
Differential age method. Every observational tests requires some 
priors, assumptions and is subjected to systematic errors. 

There is a need to develop tools which are simple and have a 
controlled systematics. The cosmological redshift drift (CRD) 
method is a test in this direction. Corasaniti et al. were the 
first to analyse dark energy models like ACDM, Chaplygin gas 
and dark energy- dark matter interaction model using the CRD 
test 1 24] . They conclude that the CRD test puts stringent con- 
straints on non-standard dark energy models. Later on, Balbi 
and Quercellini also investigated various standard and non-standard 
dark energy models ^25*] . They found the worst bound for the 
Cardassian model. Further Zhang et al. also studied Holo- 
graphic dark energy model with CRD test ||26|] and obtain a 
very tight bound on the Q^rn- The above stated work used the 
data generated by MC simulations, given by Pasquini et al. 
(2006) 

Uzan, Bernardeau and Mellier discuss the possibility in which 
the large scale structure may effect the measurements of time 
drift of the cosmological redshift 1I27I1 . Since the CRD test is 
based on the assumption that the universe is homogeneous and 
isotropic, therefore this test has been used to check the homo- 
geneity of the universe | 28]. Quartin and Amendola propose 
that CRD test can be used to distinguish between Void mod- 
els and conventional dark energy scenarios [ 291 . Recently one 
possible source of noise in the CRD test has been studied by 
Killedar and Lewis which deal with the transverse motion of 
the Lya absorbers I 3Q1 . 

Recently Liske et al. (2008) studied in detail the impact of 
next generation ELT on observing the very small redshift drift 
signal 1 6]. Using the extensive MC simulations, three sets of 
data points are generated with different strategies which include 
the measurement of the precise value of velocity shift. In our 
work. We used the recent data generated by Liske at al.(2008) 
to constrain the late time acceleration models of the universe. 
We study the models which belong to both the Einstein grav- 
ity scenario (ACDM, XCDM, scalar field potentials) and the 
modified gravity (f(R)). 

Results are summarized as follows: 



1. 



In ACDM model (Wx = -1), the;^^ minimum lies at 
fltno = 0.3. Considering /z to be a nuisance parameter, 
we marginalize over h to obtain the probability distribu- 
tion function defined as: 



Figure 4: Variation of velocity drift (cm/sec/yr) with redshift z for f(R) gravity 
model. 



4. Results and Discussions 

Till now, large number of theoretical models have been pro- 
posed to explain the observed accelerated expansion of the uni- 
verse. In order to get insight of the mechanism behind this ac- 
celeration, we require complimentary observational tools. The 



L(p) = J e-^'^^^P^'^P(h)dh 



Here P(h) is the prior probability function for h which is 
assumed to be Gaussian. We find that the best fit value of 
model parameter is independent of the choice of prior. 
We performed the;^^ statistics on flat XCDM dark energy 
model as shown in Fig. [5] The best fit values lies at 
= 0.3 and = -I. AX Icr-level the constraints are. 



30+^-^^ 



< -0.62 



5 



Here again we marginalize over h to obtain the best fit 
value of model parameters (D^, Wx). The best fit values 
are independent of the choice of prior. The redshift drift 
test give very tight constraint on the and weak bound 
on the equation of state, cOx- This is expected since the 
amplitude and slope of the cosmic velocity shift w.r.t the 
redshift is very sensitive to the ^Im and shows weak de- 
pendence on the equation of state. The other important 
feature of this analysis is that it is complementary to other 
probes such as CMB, BAO and weak lensing. 
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Figure 5: Flat dark energy model with constant cjx'. Contours in and cjx 
plane for CRD test. The best fit value lies at = 0-3 and Ux = -1. The inner 
and outer curves are at Icr and 2cr respectively. 

3 . We analyze the f(R) gravity model proposed by Starobin- 
sky with the CRD test. The variation of velocity drift 
for diff'erent values of n with redshift z is shown in Fig. 
m The expansion history of this model match exactly 
with the ACDM model for ^ = 2 and A = 2. In this 
model the Hubble parameter, H(z), becomes independent 
of the parameter n after the redshift z = 1.8 and expan- 
sion history traces exactly the ACDM model behavior at 
high redshift Jsil [23|] . Since the simulated data points 
have z > 1.9, the;^^ for this model will be the same as 
for the ACDM model. Here ACDM model means stan- 
dard cosmological model with ojx = -1, QmO = 0.3 and 
Qao = 0.7 . The eff'ective equation of state for this model 
appproaches oj = -I at the present epoch [see Fig. 8 of 
the ref.|23]]. In order to constrain this model better we 
need redshift drift data in the redshift range z< 1.9 . 

4. In Fig. O we plot the variation of x'^ with the parame- 
ter p for the exponential potential. The best fit value of 
JS = n?^J{^nf^) = 3.2 corresponds to Q^^^o = 0.26. The;^^ 
per degree of freedom is 0.67. At 3(T limit, the allowed 
range of /3 = 3.2;^-45 gives = 0.26;^-^^. Although 
the model predicts the observed value of but still this 
model is not fully in concordance with the observations; 
as the allowed range of JS at 3cr level shows that the uni- 
verse mostly remains in a decelerating phase (see Fig. O. 

5. In the inverse power law potential model, for y ^ 0, 
the energy-momentum tensor approaches towards cos- 
mological constant. The best fit values for this model are 
y = 0.025 and = 0.1 18. In Fig. |7]the contour is plotted 
in 7 - ^ plane. We would like to note here that computa- 



- exponential potential 




Figure 6: Variation ofx^ with j3 

tional eff'orts to generate the data points in the neighbor- 
hood of 7 = increases. The important feature of this 
contour is that it is complementary to the contour obtain 
by using the Dark Energy Task Force (DETF) simulated 
data sets for future experiments. This includes SNe la, 
BAO, weak Lensing and CMB (PLANCK) observations 
S. 

In the Fig. [8] variation of instantaneous equation of state 
is plotted w.r.t. redshift for the best fit values of the model 
parameters. There is almost no variation in the equation 
of state with redshift. It always stay close to 6l) = - 1 for 
the redshift range studied here. 
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k 

Figure 7: Contour iny-K plane for inverse power law potential. The inner and 
the outer contours are drawn at the Icr and 2cr level respectively. 

6. For the PNGB model, the contour is plotted in a - m 
plane as shown in the Fig. [9j The minimum lies 
at Of = 0.21 for m = 0.56. At 2cr level, a > 0.15 is 
allowed. This result is completely in concordance with 
the other cosmological observations (eg. SNe la, lensing, 
cluster observations) 1 15, 17, 18]. Again this contour is 
complementary to the contour obtain by using Markov 
Chain Monte Carlo analysis (MCMC) to generate the fu- 
ture data set. This includes future SNe la, BAO, weak 
Lensing and cosmic microwave background observations 



6 



-0.999998 r 



-0.999999 



- 1.000000 L 



2 4 6 8 

Figure 8: Variation of cj w.r.t. redshift for inverse power law potential. 
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Figure 10: Variation of u w.r.t. redshift for for PNGB potential. 



In the Fig. [TOl instantaneous equation of state is plotted 
for the best fit values. In contrast to power law potential, 
the PNGB show large variation in equation of state for 
the redshift z<2. 
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Figure 9: Contour ma - m plane for PNGB potential. The inner and the outer 
contours are drawn at the Icr and 2cr level respectively. 

It is clear that CRD test is a very simple, straightforward 
and powerful tool which probes the expansion history of 
the universe directly. In future the precision data espe- 
cially in redshift range of z < 2 shall add to the predictive 
power of CRD test. 
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